Measurement of electron wave functions and confining potentials via photoemission 
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Wave functions and electron potentials of laterally-confined surface states are determined exper- 
imentally by means of photoemission from stepped Au(lll) surfaces. Using an iterative formalism 
borrowed from x-ray diffraction, we retrieve the real-space wave functions from the Fourier trans- 
form of their momentum representations, whose absolute values in turn are directly measured by 
angle-resolved photoemission. The effective confining potential is then obtained by introducing the 
wave functions into Schrodinger's equation. 
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The electron wave functions and potentials in low- 
dimensional systems are of primary importance for tailor- 
ing electronic properties of nanostructures. The electron 
energy levels and the probability density are the phys- 
ical observables. These can be obtained in real space 
from local conductance maps via scanning tunneling mi- 
croscopy/spectroscopy (STM/STS) [Q, 0, or in reciprocal 
space by angle-resolved photoemission [|| |j . The latter 
has been used to study thin films jj| or two-dimensional 
arrays of nano-objects Q, whereas STM/STS resolves 
individual nanostructures on a surface. 

The question arises whether it is possible to directly 
derive wave functions and/or electron potentials from 
experimental data. The standard procedure consists of 
assuming a model potential, whose parameters are ob- 
tained by fitting the experiment. For instance, in core- 
level photoelectron diffraction muffin-tin potentials are 
assumed for retrieving atomic positions that fit the ex- 
periment |). The electron density of surface states 
confined between two steps has been measured by STM 
jl], |^] and has been modeled by semitransparent mirrors 
located at the step edges ||. It has been argued that 
the position of the mirror has to be moved slightly away 
from the steps to account for electron spill over . 

We introduce a direct method to determine the effec- 
tive electron potential from angle-resolved photoemission 
data, without any a priori knowledge about its nature. 
All the necessary information is contained in the momen- 
tum distribution of the photoemission intensity. This al- 
lows us to unambiguously determine the potential and 
wave functions of one-dimensional quantum well states 
on stepped Au(lll). The real-space wave functions are 
derived by Fourier transform of their momentum-space 
representations, the square of which is proportional to the 
photoelectron intensity under the conditions discussed 
below. However, the phase in momentum space is not 
measured and this leads to the well-known phase prob- 
lem in optics, shared with many other techniques, such as 
x-ray and electron diffraction |8| . Various iteration meth- 
ods have been devised to retrieve real-space objects from 
the modulus of their momentum space representation, in 



particular the oversampling method |9[). These methods 
are valid for a confined wave function, because the phase 
is obtained by repeatedly diminishing the amplitude of 
the wave function outside the confinement region. The 
one-electron potential can be obtained after dividing the 
Schrodinger equation by the wave function. 

As a test system we use vicinal noble metal surfaces 
characterized by equally-spaced, linear step arrays and 
one-dimensional surface states confined within individ- 
ual terraces of nominal size L. That is the case of 
Au(lll) vicinal surfaces with relatively wide terraces, 
such as Au(788) (L = 38 A) and Au(ll, 12,12) (L = 56 
A) shown in Fig. j§ [h}. On Au(ll,12,12), the low- 
est three quantum-well levels (N = 1 — 3) lie below the 
Fermi energy and they can be probed by photoemission, 
as shown in Figs. ||(a)-(b). The emission angle 8 is given 
relative to the terrace normal in the plane perpendicu- 
lar to the step array. The data have been taken at the 
PGM beam line in the Synchrotron Radiation Center 
(SRC) of the UW-Madison, using a hemispherical Sci- 
enta SES200 spectrometer. The energy and angular res- 
olution were 20+7 meV (photons + electrons) and 0.3°, 
respectively. The light is p-polarized and incident 60° 
off-normal. By tuning the photon energy to 60 eV we 
avoid overlap between surface umklapp replicas Jl^] . Fig. 
|2|(c) shows the angular photoemission intensity scan for 
each quantum level as obtained from the areas of the 
peaks in Fig. ||(a). The scale for the component of the 
photoelectron wave vector parallel to the terrace q x has 
been determined from the emission angle and the kinetic 
energy E as q x — sin 6 V '2m E/h. The parallel momen- 
tum scale is shifted by the reciprocal step lattice vector 
g = 27r/Z/ = 0.11 A -1 to bring it into the first Brillouin 
zone |l0| and corrected for the small parallel photon mo- 
mentum of 0.026 A" 1 . 

The photoemission intensity in reciprocal space q can 
be understood in the framework of standard photoemis- 
sion theory, assuming a one-electron description. We 
start with the photoemission intensity, which is propor- 
tional to the matrix element 

^iv(q||) oc | < q|e ik ' r A ■ V|*/v > | 2 - (1) 
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is the initial surface state wave function, |q > is 
the final electron state of momentum q = (qiutfe), qy 
is the momentum component parallel to the surface, A 
is the light polarization vector, k is the momentum of 
the photon, and N refers to the quantum number of the 
wave function. To proceed further, it will be assumed 
that the surface state wave functions can be factorizcd 
into components that are either parallel or perpendicular 
to the surface, i.e. pl| : 

*jv(r) = <f> N (R)tp{z), 

where R denotes the coordinates along the surface. If, 
in addition to the surface-state confinement in the per- 
pendicular direction z, the electron is also bound to a 2D 
region of the surface, the parallel component of the wave 
function (/>at(R) must adopt a form which depends on the 
detailed shape of the confining region and on the bound- 
ary conditions at its border. Thus, In can be written 
as 

JAr(q||)oc C\ < q|| \4>n > | 2 , 

where 

C= | <q z \e ik '*A z -^\<p> | 2 

depends very weakly on qy , so that it can be absorbed 
into a normalization factor that we set to 1. Therefore, 
there is a direct relation between photoemission intensity 
TjV-(qii) and the Fourier transform of the wave function 

<Mq||): 

<Mqn) =< qyl^iv >= / dRe- iq "- R <MR)- (2) 

Thereby we neglect multiple scattering in the final state, 
which corresponds to Fourier components q||+g|| with 
a reciprocal lattice vector gii. These affect the intensity 
distribution between different Brillouin zones but not the 
momentum distribution within one Brillouin zone (com- 
pare this to the spot profiles in LEED which are not 
affected by multiple scattering [Q). This permits ob- 
taining the wave function in real space </>jv(R) from the 
intensity in reciprocal space as 

<MR) = / ^e^-^MqiiK^i). (3) 

However, this equation still contains an unknown phase 
<5jv(q||). We have explored two methods to obtain this 
phase: (1) an iterative procedure using oversampling 
[g, g] and (2) an expansion of the wave function into a 
Fourier series combined with a least squares fit. In both 
cases, the strategy is minimizing \(f>\ outside the confine- 
ment region. 

For the iterative method, we begin with a constant 
phase <5Ar(q||) = to start the iteration with Eq. (||). 
The resulting </>jv(R) is corrected outside the confining 



region, and transformed back into qii space using Eq. 

(||). The phase of 0jv(q||) is extracted and inserted again 
into Eq. (g) to start a new iteration. The confinement 
length is obtained directly from the photoemission data 
(see below), and the noted correction consists of sub- 
tracting from the newly calculated real-space wave func- 
tion the wave function at the previous step multiplied by 
a factor 0.1. This ensures convergence to a wave func- 
tion that vanishes outside the confinement region. As a 
strong test of convergence of this iterative technique, we 
have used the phase of the N — 2 state to start the iter- 
ation of the N = 1 state. The first iteration step leads to 
a wave function that resembles that of the N = 2 state, 
a clear indication of the importance of the phase |g] , and 
convergence to the true N = 1 state is achieved after 
several hundred iterations. 

This iterative method leads to results that are in agree- 
ment with those derived from an expansion of the wave 
function into sine functions that vanish at the edges of 
the confinement region. Moreover, the length of the re- 
gion where the wave function takes non-negligible values 
is not sensitive to the input value for the length of the 
confinement region. 

The size of the confining region is obtained from the 
Fourier transform of the experimental intensity distribu- 
tion, which can be written as the self-convolution of the 
real-space wave function, 

/ (2^ eiq "' R/Ar(q 'i ) = / dR '<MR')<^(R'-R)-(4) 

This convolution takes non-zero values in an area twice 
the size of the confining region. For the data of Fig. || in 
particular, one obtains curves that are basically confined 
within a region of « 60 Ain diameter (not shown), in 
agreement with the nominal terrace width in Au(ll 12 
12), L = 56 A. 

The natural normalization for both the two- 
dimensional wave function 4>n and the measured inten- 
sity is provided by Eq. ([|) if one sets R = 0: 




Zjv(qn) =< <M0w >= L 



Applying this procedure to the angular scans of Fig. ||(c), 
we have obtained the surface state wave functions for 
quantum- well levels in Au(ll 12 12) that are shown in 
Fig. [I]. All three wave functions are confined to a region 
whose width matches the terrace width L. Thus the wave 
functions exhibit a clear terrace confinement that has 
not been assumed by our reconstruction procedure, but 
rather it has emerged from the information contained in 
the photoemission data. 

The actual effective potential of the terrace can be re- 
trieved from the Schrodingcr equation as 

nR)-^ = ^^v 2 MR), (5) 



3 



where m* — 0.26m is the effective mass of the electron 
in the initial state of energy Ejq , taken from the disper- 
sion of the surface state on flat Au(lll). The potential 
becomes independent of the quantum number N if many- 
electron effects are absent, as we assumed initially. 

The experimental wave functions of Fig. [I] have been 
introduced into Eq. (|B|) and the resulting electron poten- 
tial has been represented in Fig. [|. In order to compare 
the potential derived from the quantum states with var- 
ious N we add on the left side of Eq. (||) the respec- 
tive experimental values E^ p shown in Fig. |l| (horizon- 
tal bars). Excellent mutual agreement is obtained for the 
shape of the potential derived independently from each 
of the three wave functions N = 1 — 3, proving the va- 
lidity of our method. The potential exhibits a smooth 
central region and sharp boundaries that force electron 
confinement. It must be stressed that the validity of Eq. 

is limited to regions where the wave function is not 
too small, and therefore, the asymptotic limit of the step 
barrier potential cannot be determined. 

There is a slight uphill/downhill asymmetry of <f> and 
V, which is related to an asymmetry of \4>\ for ±q . This 
effect shows up for the highest level (JV = 3), which 
is the least confined. Such an asymmetry explains the 
asymmetric reflectivity observed in STM Q where the 
downhill step reflects more strongly than the uphill. The 
asymmetry increases for higher-lying empty states which 
start spilling over the weaker potential barrier on the up- 
hill side. 

The same iterative procedure has been applied to the 
one-dimensional quantum well of Au(788), measured in a 
different system Q . In this case the terraces are 38 Awide 
and only the first two quantum levels are occupied. We 



also obtain a good consistency in the electron potential 
between the two levels. The resulting average is displayed 
in Fig. |l[ and compared with the average potential of 
Au(ll,12,12). In contrast to the latter case, the width of 
the potential well in Au(788) is narrower than the terrace 
size, as also observed in thin Pb films 0. 

In summary, we have introduced a simple procedure 
for directly obtaining wave functions and effective poten- 
tials of confined electronic states from the momentum 
distribution of the photoemission intensity. This pro- 
cedure, which relies on photoemission theory within the 
one-electron approach, has been successfully tested in the 
case of laterally confined surface states on Au(lll) vic- 
inal surfaces. The method is of general applicability for 
discrete states because they are spatially confined and 
allow an iterative determination of the phase by forcing 
0=0 outside the confinement region. It can be applied to 
any type of surface nanostructure confined in at least one 
dimension, such as arrays of quantum wires and quan- 
tum dots. While STM provides direct images of the low 
spatial frequencies, photoemission determines the high 
spatial frequencies. It would be interesting to pursue the 
reverse transformation from the real space \<f>\ in STM to 
<f> in q space using the technique described here. A con- 
nection between \<fi\ and the E(q) band dispersion has 
been made already [l3| 
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FIG. 1: Vicinal Au(lll) surfaces with quantum well states 
confined by step edges. The energy levels are directly mea- 
sured by angle- resolved photoemission. The confining poten- 
tial V and the wave functions of states N=l ,2,3 are obtained 
from the momentum distribution of the photoemission inten- 
sity. 
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FIG. 2: (a) Photoemission spectra from Au(ll,12,12) show- 
ing one-dimensional surface states, (b) Second energy deriva- 
tive of the intensity showing the presence of three quantum 
levels E^ F . (c) Peak area of each quantum level obtained by 
line fitting of individual spectra in (a). The lines are spline- 
fits to the data points. 




FIG. 3: Confining potentials V obtained from the experimen- 
tal photoemission intensities of Fig. ^ for the JV=1,2,3 states 
in Au(ll,12,12) using the oversampling method. 



